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ABSTRACT 

We present a 5' x 5' integrated intensity map of ^^CO (J = 3 — 2) emission from the p Ophiuchi 
cloud core that traces low-luminosity outflow emission from two protostars: Elias 29 and, most 
likely, LFAM 26. The morphology of the outflow from Elias 29 is bipolar and has a curved axis 
that traces the S-shaped symmetry seen in H2 emission. The outflow from LFAM 26 is a new 
detection and oriented in the east/ west direction near the plane of the sky with most of the 
blue-shifted emission being absorbed by intervening clouds. The outflow axis of this object also 
appears to intersect a knot of H2 emission previously attributed to Elias 29. LFAM 26 is a low 
luminosity source (iboi = 0.06 L©) which, in combination with the observed outflow, makes it 
a candidate Very Low Luminosity Object (VeLLO). We derive lower limits to the gas column 
densities and energetics for both outflows. The mechanical luminosities for Elias 29 and LFAM 26 
are 6.4 and 10.3 xlO^'^ i©, respectively. 
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1. Introduction 

We present ^^CO (J = 3 — 2) observations towards the (36 L©) Class I protostar Elias 29 (a = 
16'>27'"09".5, 5 = -24°37'18" [ J2000.01) locat ed in the p Ophiuchi star forming region. The p-Oph cloud 
complex is nearby {d ~ 120 pc, lLoinardll2006l ) and has a diverse (proto-)stellar environment that makes it 
a superb laboratory for studying the impact of molecular outflows on the surrounding environment. The 
region around Elias 29 is a particularly well-studied component of the cloud complex. Two foreground clouds 
at vlsr = 2.7 km s~^ and vlsr = 3.8 km s~^ account for the presen ce of cool ice along the hue of sight 
and absorb much of the blue-shifted outflow emission in this region (jBoogert et al.l l2002f ) . Furthermore, 
Boogert et al. ( 20021 ) i dentified a ridge o f den se material at vlsr = 5-0 kms~^ within which both Elias 29 



and LFAM 26 reside. 



Ceccarelli et al.l (j2002r ) obtained spectra of the ^^CO ( J = 6 — 5) transition with 
the JCMT and the J > 15 CO levels with the ISO SWS towards Elias 29. These authors found evidence 
for a super-heated su rface disk layer with a temperature and mass similar to that found in Herbig AeBe 
stars (jThi et al.ll200l[ ). suggesting that Elias 29 may be a deeply embedded Herbi g AeBe star or a transition 
object between the standard Class I and Class II types. iBontemps et al.l (|l996f ) reported the discovery of 
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an outflow around Elias 29 and estimated its moment um rate. A larger map of the region was made in the 
12C0 (J = 2 - 1) transition bv lSekimoto et all |l997h . 

In this paper, we present ^^CO {J — 3 ~ 2) observations that further expand the size of the mapped 
region and reaffirm the bipolar nature of the outffow, but find that it is more extended than previously 
thought. Furthermore, the map has sufficient size to identi fy a previously und iscovered outfiow whose source 
appears to be LFAM 26. First discovered as a radio source (jLeous et al.lll991l ). LFAM 26 w as later identified 
as a low luminosity (Lboi — 0.064 X©) Class I protostar based on ISOCAM photometry (jBontemps et al 



20011 ). The extreme faintness of this source, coupled with the detection of molecular outflows presented in 



this paper, make it a candidate VeLLO. However, detailed modeling of the infrared SED will be necessary 
to determine if this object meets the primary criterion of Lint < 0.1 Lq. 



Observations 



We made a 5' x 5' ^^CO ( J = 3-2) On-The-Fly map centered spectrally at 345.795975 GHz and spatially 
on Elias 29 using the 10m Heinrich Hertz Telescope (HHT) on Mt. Graham, Arizona during April 2001 and 
January 2006. A 256 channel filterbank spectrometer with 250 kHz resolution was used as the backend in 
2001, and two 2048 channel Acousto-Optical Spectrometers ( AOS) with 475 kHz resolution were used in 2006. 
The off-position for both maps was located at a clean offset ( Kulesa et al. 20051 ) of Aa = 2.5°, = —0.83° 
from Elias 29. We obtained two rows of data per OFF, translating to ~1 OFF observation per minute for 
a 10" per second scan rate. Each row is offset in declination by 10". At 345 GHz, the half-power telescope 
beamsize and main beam efficiency were measured to be 22" and 0.7, respectively. Pointing accuracy of a 
few arcseconds was maintained by observing planets approximately once per hour. The OTF data grid was 
convolved with an 11" beam to obtain a Nyquist-sampled map. System temperatures during the observing 
run were 500-1500 K. Line stre ngths were calibra ted using hot and cold loads in conjunction with sky 
measurements at the off position (jStark et al.ll200l[ ). Calibration uncertainty is estimated to be < 20%. For 
both datasets, the main beam temperature (T^b = T\/r]mb) showed a typical rms noise level of 0.7 K . The 
2001 data are smoothed in velocity space to match the 0.4 km s~^ bins of the 2006 data. The final datacube 
used in the analysis is an rms- weighted average of the two datasets and has a typical rms noise temperature 
of T^fc = 0.4 K. 



Results 



3.1. Outflow Morphology 



In Figure [1] we show spectra located at the driving sources as well as at the approximate peak of the 
integrated emission intensity in the blue and red lobes of both outflows. The vertical lines demarcate the 
position in velocity space of intervening absorption clouds along the line of sight at f r,,SB = 2 .7 kms~^ and 
vlsb. = 3.8 km s"-"^, as well as a ridge of material at vlsr = 5 km s"-"^ ( Boogert et al. 20021 ). Each of the 
spectra show signs of signiflcant absorption due to these intervening clouds. 

Figure [2] shows the integrated intensity map of the blue-shifted (solid contours, vlsr — —4 to 1 km s^^) 
and red-shifted (dashed contours, vlsr = 7 to 12 km s^^) outflow emission. Th e velocity rang e s hav e 
been chosen to avoid the absorption features from the intervening clouds detected bv lBoogert et al.l (j2002l ). 
Crosses show the location of Elias 29 and LFAM 26, plus signs indicate the position of infrared sources 
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(|Wilking fc Ladalll983h . s quares represe nt starless cores (IMotte et al.l 119981) . and the diarnonds show H2 
emission knots observed by Gomez et aL (feature G2 2003 ) and by Ybarra et al j (features l-3b 2006 ). Ehas 29 
is associated with emission tow ards the north (blue-shifted) and south (red- shifted) , consistent with results 



from previous outflow surveys (jBontemps et al 



outflow with a mostly pole-on inclination (jHanner et al.lll995l ). 



1996 



Sekimoto et alJ 119971 ) that have suggested a bipolar 



The large scale outfl ow morphology stro ngly resembles the S-shaped symmetry observed in H2 emission 
near the driving source (| Ybarra et alJ 120061 ). with the blue lobe bending to the northeast and the red lobe 
bending to the southwest. The solid curve through the peaks of the Elias 29 blue and red lobes represents 
our best estimate of the outflow a xis and is weighted by eye towards the five H2 emission knots labelled 
as features 1, 2a, 2b, 3 a, and 3b (Ybarra et al. 200(k i in the main panel of Figure [21 The bottom right 



, , ahM—^^^MMMr^^^''^ ' V i: O t=J --0 

inset is reproduced from I Ybarra et al.l ( 20061 ) and shows the structure of the region as seen in H2 em ission. 



Regions 3a and 36 appear to trace the peak of the red-shifted outflow emission. ICeccarelli et al ] (I2OO2) 
mapped CO ( J = 6 — 5) line emission over a 35" x 35" field centered on Elias 29 with an effective beamsize 
of 12". Contour maps of their observed wing emissio n show evidence in support of an east-west outflow 
axis, consistent with the S-shaped symmetry found by Ybarra et al. ( 2006h and seen in the central regions 
of our map. Higher spatial resolution observations will be necessary to explore the detailed morphology of 
the central region of outflow emission from Elias 29 and provide a more precise comparison with the known 
H2 regions. 

In addition to the well-studied o utflow from Elias 29, we report a new out flow whose source appears 



to be LFAM 26, a Class I protostar (jOreene et al.lll994l : iBontemps et al.l 120011 ) located ~1' northwest of 
Elias 29. The morphology of this outflow — with lobes to the east and west both showing signiflcant red- 
shifted emission — suggests that its axis lies near the plane of the sky, as shown by the solid line in Figure [H 
The relatively low level of blue-shifted emission from this outflow is mo st likely due to the presence of the 
foreground absorbing clouds at vlsr = 2.7 kms~^ and vlsr = 3.8 kms~^ (jBoogert et al I2OO2I) . In particular, 
the Vlsr — 2.7 kms^^ cloud is roughly spatially coincident with the eastern lobe of the LFAM 26 outflow. 
H2 em ission detected within the contours of the integrated emission intensity to the east (jYbarra fc Barsonv 
20041 ) lends further support to the notion of LFAM 26 as the driving source. The outfl ow axis as drawn lie s 
nearly a full beamwidth away from [GSWC2003] 2, a knot of II2 emission detected bv lGomez et al.l (|2003f ). 
These authors attributed the H2 knot to the Elias 29 outflow, but in our CO ( J = 3 — 2) wing map, the 
western lobe of the LFAM 26 outflow appears to he between LFAM 26 and [GSWC2003] 2, suggesting that 
LFAM 26 is a more likely candidate than Elias 29 to be the powering source of the II2 knot. 

It must be kept in mind that this is a particularly dense region of p-Oph which makes the association 
of outflows with their driving sources particularly problematic. For this reason, additional evidence linking 
LFAM 26 and the observed outflows is desirable. This evidence comes in the form of nebulous emission 
seen in deep NIR imaging of LFA M 26 that has a dist inct hourglass shape — shown in the top right inset 
of Figure (taken from Fig. 4 in iDuchene et al.l 120041 ) — with the nebulosity extending a few arcseconds 
away on either side of the central source. The orientation of the hourglass roughly coincides with the ~100° 
outflow position angle we estimated above, serving as further evidence that LFAM 26 is the powering source 
of the observed outflow. The hourglass shape around LFAM 26 is consistent with the presence of an edge-on 
accretion disk. Whether or not LFAM 26 qualifles as a VcLLO must await detailed modeling of its IR SED, 
a task which is beyond the scope of this letter. 

It should be noted that the integrated emission intensity peaks do not appear to lie directly on the 
outflow axis as shown in Figure [2l Due to the masking effect of foreground clouds, the redshifted emission 
components alone may not accurately portray the east-west spatial centroids of the outflow. Nevertheless, let 
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us consider the possibility of alternative driving source(s). First, the object closest to LFAM 2 6 is E-MM3 
located ~15" east and ~30" south of LFAM 26. This object was first labelled as a starless core (jMotte et al 



19981 ). but subsequent VLT H and Ks o bservations showed a bipolar reflection nebulosity surrounding it 



in much the same manner as LFAM 26 (|Brandner et al.l 12000) . However, the orientation of this reflection 
nebula is closer to north/south than east /west and is therefore an unhkely candidate for the driving source. 
Another possibility is that each of the outflows is driven by one of the other embedded objects either in this 
region (e.g. E-MM4 or E-MM5) or beyond. Given the highly obscured, complex structure of the region, 
this scenario is plausible. Finally, it is conceivable that Elias 29 is a binary system, with the other object 
producing th e observed outflow . However, numerous se arches using a wide v ariety of techniques from lunar 
occultations (jSimon et al.l 119951 ) to deep NIR imaging (jPuchene et al.l 120041 ) have found no evidence for a 
companion. 



3.2. Outflow Energetics 



Physical parameters of each outflow are summarized in Table [T] From the integrated intensity map, we 
estimate the maximum radius {Rmax) and area {A) of each outflow lobe centered on the position of peak 
flux. Here the characteristic velocity of each outflow, Vchar, is the maximum outflow velocity detected in 
the associated ^^CO (J = 3 — 2) spectra, ~ 8kms^^. The dynamical timescale is t^i — Rmax/^char- We 
assume the wings of the ^^CO (J = 3 — 2) emission are optically thin, so the values in Table [1] represent lower 
limits to the actual values. As mentioned in Section 13.11 however, there is evidence of strong absorption 
of the blue-shifted outflow emission from LFAM 26 due to the presence of foreground absorption clouds. 
Furthermore, the measured radial velocity underestimates the characteristic velocity by a factor of l/sin{i), 
where i is the inclination angle. If the outflow axis of LFAM 26 is near the plane of the sky, then i approaches 
0° and the correction factor can increase significantly. 

From the peak main beam temperature in the ^^CO ( J = 3 — 2) spectrum, we estimate T^x — 25 K. 
While Tex is not constant across the map, the variations are small and do not significantly affect the results. 
For instance, if we assume a T^x that is smaller by a factor of 2, we find average c olumn densiti e s that 
are roughly a factor of 2 larger. The total CO column density, Nt, is calculated from IXulesa et aP ( 20051 ) 
under the assumption of LTE. Assuming dark cloud abundances of [H2]/[C0]= 1 x 10^, we estimate the 
mass of each outflow lobe to be M np = ([H2/[CO])AfT^Aigm-H2 1 where /ig is the mean atomic weight of 
the gas, 1.36 (jSekimoto et al.lll997l ). The mass outflow rate is Mqf = Mop/td, the outflow momentum 
is Pop = MopVchar and the outflow force is Fqp — Pop /id- The mechanical luminosity is estimated as 
Lm = \Mopv'qp /td- The outflow to infall conversion efficiency is e = Mqp / Maccr- 

Sekimoto et al.l ( 1997 ) and Bontemps et al. ( 19961 ) have obtained integrated intensity ^^CO (J = 2 — 1) 
emission maps centered on Elias 29. The ^^CO (J = 3 — 2) map shown here represents a significant 
enlargement of the surveyed area. Inside t he mapped region, we find a larger radius for the blue lobe of 
Elias 29 — 0.058 pc compared to 0.024 pc in lSekimoto et al" ( 1997 ). We also find an average column density 
over the Elias 29 outflow that is a factor of ~2 — 3 lower than observed in ^^CO ( J = 2 — 1). The difference in 
column density estimates may be due to the ^^CO (3 — 2) wing emission being optically thick. Our optically 
thin assumption would then underestimate the true column density by a factor of t/(1 — e""^). Due to the 
extended size o f the o utflow region, the ^^CO ( J = 3 — 2) derived outflow mass is comparable to that of 
Sekimoto et al.l (jl997[ ). The outflow energetics are summarized in Table [TJ The low values of the energetic 
parameters ass ociated with Elias 29 s uggest it is near the end of its outflow/ accretion phase, consistent with 
the findings of ICeccarelli et al.l (|2002r ). 
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4. Summary 

We present a 25 square arcminute map of ^^CO ( J = 3 — 2) outflow emission centered on Elias 29. In 
the map. wc see evidence suggesting the Elias 29 outflow is bipolar and has a curved flow that traces the 
same structure seen in H2 line emission. We report the first detection of a molecular outflow from LFAM 26, 
an object whose low luminosity makes it a candidate VeLLO. While the high density of protostellar objects 
in this region produces a complicated outflow structure that makes associating outflows with their driving 
sources somewhat problematic, we conclude LFAM 26 is the most likely driving source. Both the molecular 
outflow and the NIR scattered light cavity morphology suggests the associated outflow axis lies near the plane 
of the sky, in which case we are viewing the protostar through an edge-on accretion disk. We calculate lower 
limits on the mechanical luminosity of the outflows from Elias 29 and LFAM 26 to be 6.4 and 10.3 x 10~^ Lq, 
respectively. The low value of for Elias 29 suggests this object is near the end of its outflow/accretion 
phase of evolution. 

We wish to acknowledge the comments of the referee which were very helpful in refining the paper. 
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Table 1. Parameters of Detected Sources 



Parameter 


Elias 29 


LFAM 26 


■^max (P^) 


0.058 


0.064 


td (103 yr) 


7.1 


7.8 


<Nco> (10^^ cin-2) 


3.6 


3.3 


^OF (sq. arcsec) 


5100 


6600 


MoF (IO-^Mq) 


8.7 


15.2 


MoF (IO-6M0 yr-i) 


1.2 


1.9 


e {MoF/Maccr) 


0.15 


0.25 


PoF (lO-3M0knis-i) 


69 


121 


FoF (lO-'^Mokms-i ) 


9.8 


15.6 


Lm (10-3 Lq) 


6.4 


10.3 
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n 
E 




-5 5 10 

Vlsr (km s-1) 



15 -5 5 10 

Vlsr (km s-') 



15 



Fig. 1. — Spectra of Elias 29 and LFAM 26 in addition to their respective outflows. The spatial offsets from 
the driving sources are given below each spectrum (north and east are positive). LFAM 26 is offset from 
Elias 29 by (-60", 44"). Dotted lines indicate the velocity of intervening absorption clo uds at vlsb^ ~ ^■'^ ^^ '^ 
VLSR = 3.8 kms~^ as well as an absorption ridge at vlsr = 5.0 kms~^ detected bv lBoogert et al.l (|2002l ). 
The spectra have been shifted in intensity by 15 K. 
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Fig. 2. — Integrated intensity map of the blue-shifted (sohd contours, vlsr = —4 to 1 km s~l) and the 
red-shifted (dashed contours, vlsb. = 7 to 12 kms^^) outflow emission centered on Ehas 29. Lowest contour 
shown represents 3cr, steps are in units of 2a. cfby^ = cr,.^ = 1.4Kkms~^ . Crosses show the lo cation of 
Elias 29 and LFAM 26, plus signs indicat e the position of infrared sources dWilking fc Ladalll983h. s quares 
repr esent starless c ores (Motte et al. 19981 ). and diamonds show H2 emission knots from lComez et al. ( 20031 ) 
and lYbarra et al.l (j2006| ). A solid line intersecting LFAM 26 is drawn to indicate the ~100° position angle 
of the LF AM 26 outflow and is reproduced in the top right inset, which shows continuum H2 emission from 
LFAM 26 lOuchene et al.ll2004l ). This hue is nearl y a beamwidth away from [GSWC2003] 2 (designated G 2 
in the figure), an H2 emission feature detected bv lGomez et al.l (|2003l ). A solid curve through the peaks of 
the Elias 29 outflow lo bes indicates the S- shaped symmetry of the outflow and is weighted by eye to trace 
the H2 emission from (jYbarra et al.l 120061 ). The bottom right inset shows H2 emission and a lOcr contour 
of continuum emiss ion detected in the narrowband H2 filter used for the observations. The physical scale is 
shown in the inset |Ybarra et al.ll2006h . The 22" FWHM beam is shown m the top left corner. 



